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^ ' We propose new effective inter-nucleon forces with a finite-range three-body operator. The 

^ I proposed forces are suitable for describing the nuclear structure properties over a wide mass number 

region, including the saturation point of nuclear matter. The forces are applied to microscopic 
calculations of ^ = A'^ (A < 40) nuclei and O isotopes with a method of antisymmetrized molecular 
' dynamics. We present the characteristics of the forces and discuss the importance of the finite-range 

(N , three-body term. 
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I. INTRODUCTION 
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' Experimental information concerning nuclear structure has been rapidly increasing in the region of unstable nuclei 
, as well as stable nuclei. Together with the increase of experimental data, microscopic calculations of nuclear structure 
have developed to describe those structure properties. For a systematic study of various types of structures covering 
a wide region of the nuclear chart, the key issue in microscopic approaches is an extension of the model space. 
Corresponding to the advancement of microscopic calculations, one must improve the effective forces so as to be 
^ I suitable for extended model space. 
' One of the basic properties of nuclei is saturation of the binding energy and density. The saturation property 

implies that a nucleus can be divided into fractional nuclei, and that such nuclei can fuse into a unified nucleus with 
^ . only a small amount of energy. This "easiness of assembling and disassembling nucleons" [1] is one of the fundamental 
^ ' features of the nuclear many-body system. In nuclear-structure problems, this dynamics appears as a variety of 
^ , structures such as clustering of nucleons and mean-field formation. Indeed, it is known that there exist many cluster 
states not only in stable nuclei, but also in unstable nuclei. This means that the basic dynamics of the nuclear 
structure contains two very different fundamentals. One is the formation of clusters, and the other is the formation 
of a mean field. With the traditional microscopic approaches, it was difficult to describe both of the features in one 
" " " framework because of the limitation of the model space. Typical mean-field approaches, such as Ilartree-Fock(HF) 
methods, are applicable to studying the mean-field aspect in the heavy nuclear region, while they are not suitable for 
studying cluster structure in light nuclei. On the other hand, cluster features in light stable nuclei were studied by 
cluster models. Recently, hybrid models, which can describe both the cluster aspect and the mean-field aspect, were 
developed and applied to structure studies of unstable nuclei as well as stable nuclei. One of the powerful approaches 
is a method of antisymmetrized molecular dynamics(AMD) [2,3], which is able to describe two different dynamics, i.e. 
the clustering and the mean-field formation due to its ab initio nature. In the study of unstable nuclei in the p-shell 
and sd-shell regions, performed with the AMD [1-3] method, it was found that both the cluster features and mean- 
field are essential in a systematic study of the ground and excited states of unstable nuclei as well as stable nuclei. 
As for an extended model based on three-dimensional HF calculations [4,5], a variation after the parity-projection 
method [6] might describe a variety of structures in light nuclei. In a systematic study with such extended models, 
the quality of the adopted effective interactions becomes an increasingly serious problem, because the structure is 
based on balance and competition between two different dynamics, like cluster formation and mean-field formation. 
This means that some theoretical results are sensitive to the adopted effective force, namely the Hamiltonian, which 
governs the system. In other words, the effective force must be appropriate to describe the characteristics of various 
states which may appear in the extended model space. 

In microscopic calculations of nuclear structure, one often uses the effective forces to solve nuclear many-body 
problems. In most of the phenomenological effective inter-nucleon forces, the dominant part is expressed by central 
forces. With the Minnesota force [7] and the Volkov force [8], which consist of two-body central forces, the properties 
of very light nuclei are well reproduced. However, it is impossible to describe the saturation property of nuclear 
matter and heavy nuclei (mass number A > 10) with such two-body central forces. For satisfying the saturation 
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properties, phenomenological density-dependent repulsive terms arc usually imported in addition to the central two- 
body terms. Such forces with a density dependence are Gogny forces [9], Skyrme-type forces [11,10] and Modified 
Volkov forces [12], for example. These forces contain zero-range density-dependent operators or zero-range three-body 
operators. HF calculations with Gogny forces and Skyrme forces have succeeded in describing the bulk properties of 
the low-lying states of hc;avy nuclei. These forces, however, fail to reproduce the properties of very light nuclei. Tohsaki 
introduced repulsive terms with finite-range three-body operators, and proposed new effective forces responsible for the 
microscopic a-cluster model, which succeeded to describe the overall properties over a wide range of masses from a to 
nuclear matter [13]. Namely, the a-a scattering behavior, the size and energy of double closed shell nuclei(Q;, '^^O and 
'^''Ca) and the matter saturation property were successfully reproduced with Tohsaki's forces. However, unfortunately, 
the forces contain too complicated three-body operators to be employed in practical microscopic calculations of general 
nuclei. 

Our aim in this paper is to propose new effective inter-nucleon forces which are applicable to practical calculations of 
general nuclei with the AMD method. In order to save computational costs, we introduce a simple three-body operator 
expressed by a one-range Gaussian. The new forces are applied to microscopic calculations oi Z = N nuclei up to *°Ca 
and neutron-rich O isotopes with the AMD method. We explain the characteristics of the proposed forces and prove 
the success in describing the overall properties from the a particle to nuclear matter, by demonstrating the results 
of the radii and energy of finite nuclei, the elastic a-a scattering phase shift, and the saturation property of nuclear 
matter. We point out the significance of finite-range three-body terms in the effective forces from a phenomenological 
point of view, and discuss the role of the repulsive three-body terms in the connection with the tensor force, which is 
the origin of matter saturation. 

This paper is organized as follows. The new effective forces with a finite-range three-body term are proposed in II. 
We describe the characteristics of the effective forces by demonstrating the calculated results with the new forces in 
III. In IV, the significant role of the finite-range three-body term is discussed. Finally, in V we summarize the present 
work. 



II. EFFECTIVE CENTRAL FORCES 

As mentioned above, for most of the phenomenological effective forces in microscopic models, the dominant part is 
the central force. Concerning the effective central forces, it is regarded as a tensor force, a hard core and other terms 
in the bare interactions are renormalized because it is not easy to explicitly treat such terms in microscopic models 
for a nuclear many-body system. By categorizing the effective central forces into three types, we briefly review the 
general tendency of the typical effective forces to reproduce the structure properties, such as the saturation behaviour 
of nuclei and nuclear matter and the a+a inter-cluster potential. 

(1) Regarding finite-range two-body forces with no density-dependent terms, there exist such effective forces as 
Minnesota [7] and Volkov [8] forces which can reproduce the a-^a phase shift and the size and binding energy of an a 
particle and light p-shell nuclei. However it is impossible to represent the saturation of matter and heavy nuclei with 
such the two-body central forces. 

(2) Gogny [9] forces, Skyrme-type [10,11] forces and Modified Volkov [12] forces, which can guarantee nuclear 
matter saturation, contain zero-range density-dependent terms or zero-range three-body terms. Gogny forces and 
Skyrme-type forces successfully explain the size and the binding energy of the overall heavy nuc;lei, and also matter 
saturation. In spite of the success of nuclear matter saturation, all of these effective forces fail to reproduce the size 
of the a particle and the a-a scattering behaviour. Namely, all of these forces give larger a-particle sizes and a less 
a-a inter-cluster potential than the experimental ones. Therefore, the applicability of such forces to very light nuclei 
is not assured. 

(3) We know only one kind of the effective force which can simultaneously guarantee the size of a particle, the 
a-a. phase shift, and the saturation properties of nuclei and nuclear matter. That is the third type of effective force 
containing finite-range three-body operators, proposed by Tohsaki [13], namely, Tohsaki's Fl and F2 forces. It was 
found that these effective forces can successfully reproduce the a-a phase shift, the sizes and binding energy of 
double-closed nuclei up to ''°Ca and also matter saturation point. 

In the present paper, we propose new effective forces which can be classified into the last category (3), forces with 
finite-range three-body operators. To reproduce the nuclear matter saturation properties, the effective central forces 
must contain density-dependent terms. It is considered that the origin of the density-dependent terms is suppression 
of the tensor force of bare nucleon-nucleon interactions in the nuclear matter. In other words, the density-dependent 
terms arc considered to phenomcnologically simulate the medium effect that tensor force is suppressed in the nuclear 
matter because of the Pauli principal. From this point of view, it seems to be reasonable that the effective central 
forces contain finite-range three-body terms because the original tensor force is the middle-range one. Moreover, the 
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TABLE I. Parameter sets of the effective central force(F3B force) in Eq. 1. 



vf^(MeV) vi^\MeV) a^^\im) a^^^(fm) 


vf)(MeV) (fin) 


v(3)(MeV) a'-^^fm") 




m 


b 


h 


-198.34 300.86 1.2 0.7 


22.5 0.9 


600 1.25 


case(l) 
case(2) 


0.421 
0.193 


0.1 
-0.185 


0.085 
0.37 



type (2) forces with zero-range density-dependent terms have a fatal difficulty in simultaneously describing the density 
and the energy of an a particle and nuclear matter. This failure in the quantitative reproduction of the saturation 
properties with type (2) forces is proved later in IV. 

We propose effective central forces, Vcentrai, expressed by two-body and three-body operators as follows: 



Vcentral — + 



^ijk ' 



i<j<k 



where 



= (1 - m + - hPr - mP^Pr) [v^^^ oxp [-(^)'] } + 4'^ exp 

= [__^ _ r .)2 + - r,)^ + (r^ - r,)^}] . 
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Here we adopt a simple operator for the finite-range three-body terms. Compared with Tohsaki's forces expressed by 
three-range three-body operators, this one-range three-body operator with a symmetric form in the present forces is 
suitable for practical microscopic calculations of general nuclei with the AMD method. 

We propose a set of parameters for the ranges and strengths of the two-body and three-body terms. The dominant 
components of the two-body terms are the first term in Eq.2. We adopt the same form as the Volkov force for this 
term, which consists of a medium-range attractive Gaussian and a short-range repulsive Gaussian. Wigner, Bertlett, 
Heisenberg and Majorana terms are chosen to be same parameters as the common multiplier of these two range 
Gaussians. The basic idea of the derivation of the parameters is that the nuclear structure properties over a wide 
mass number region should be described with a two-range two-body force and an one-range three-body force without 
mass-number depending parameters. We adopt the same range parameters for the two-range Gaussian of the two- 
body force as the Gogny force(ai^^ = 1.2(fm) and = 0.7(fm)). With longer range parameters, we failed to obtain 
sufficient reproduction of the nuclear structure. Although these range parameters for the two-body terms are shorter 
than the Volkov and MVl forces, the range of the inter-cluster potential is not so different from those given by the 
MVl forces, as shown in later, due to the effect of the repulsive finite-range three-body term in the present force. In 
order to explain the nucleon-nucleon scattering phase shift, the form with a two-range two-body and an one-range 
three-body is too simple because of the restriction of the simplified three-body operator. We add a small amount of a 
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repulsive two-body term given by the second term of the two-body force in eq.2, v\ exp 
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and tune 



again the strength of this term and the three-body term. In table I, the proposed parameter sets of the present forces 
are presented. We call the present parametrization as the FSB forces. Two sets case(l) and case(2) of parameters of 
m, b and h for Wigner, Bertlett, Heisenberg and Majorana terms are listed in table I. Comparing the case(l) and 
case(2) parameters, the triplet and singlet even components arc same but the odd components different between two 
sets (1) and (2). As a result, the case(l) and case(2) interactions give the same results for double closed Z = N nuclei 
written by HO wave functions, the a-a interaction, and total energy of symmetric nuclear matter, however, they give 
different results for Z ^ N nuclei. 

The present parameters m for the Majorana term shown in table I are determined so as to reproduce the binding 
energy oi Z = N nuclei in mass region A = 4 ^ 40 with the present AMD-I-GCM calculations. The Majorana 
parameter m is considered to be an adjusting parameter in fine tuning to fit the experimental data. It means that 
one can modify the value m, which is suitable to his model wave function. 

In practical AMD calculations of Z = N nuclei and O isotopes, we also employ the two-body spin-orbit force in 
addition to these central FSB forces. 
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FIG. 1. Calculated phase shift of elastic a-a scattering with the FSB force, the MVl case l(m = 0.60), and the Volkov No. 2 
(m = 0.60 and m = 0.65). The width paramctcr(i^) for an a particle in an RGM calculation was chosen to bo 0.257, 0.209, 
0.267 (fm~'^) for the FSB, MVl, Volkov No. 2 forces, respectively. The points indicate the experimental data taken from [16]. 



III. CHARACTERISTICS OF THE EFFECTIVE FORCES 

In this section, we exhibit the characteristics of the proposed effective forces concerning the saturation properties 
and a-a scattering behavior compared with other types of typical phenomenological effective forces. 

A. a-a interaction 

The a-a inter-cluster interaction is important for studying the cluster structure of light nuclei. The scattering phase 
shift is a useful physical quantity to inspect the features of the inter-cluster interaction. The phase shift of the elastic 
a-a. scattering is calculated by the resonating group mcthod(RGM), where an a particle is simply written using a 
(Os)'' harmonic oscillator wave function [14]. The width parameter of a clusters is chosen to be the optimum value 
for an isolate a-particle for each effective force. 

In Fig. 1, the calculated a-a phase-shift values are compared with the experimental data. It has been known that 
the experimental phase-shift data are well reproduced by the Minnesota force [7] and the Volkov No. 2 force(m = 0.60) 
[15]. The theoretical values calculated with the new force(F3B) agree reasonably well with the experimental data. 
On the other hand, the MVl casc(l) forcc(m = 0.60) underestimates the experimental data, which means that the 
force provides a less-attractive a-a interaction. It was also foimd that the SII force(a Skyrme-type force) gives a 
less-attractive inter-a interaction, and therefore, fails to reproduce the a-a phase shift [13]. 

Figure 2 presents the energy and potential as a function of inter-cluster distance in the a-a system. The a-a system 
is represented by the Brink-type wave function. Namely, two a clusters are located at points, {x,y,z) = {0,0, d/2) 
and {x,y,z) — (0,0, —d/2), where d, indicates the inter-cluster distance. The width parameter of an a is fixed to be 
v = 0.25 fm~^. All of the energy curves have minimum points at around d ~ 3 fm. The energy and potential for F3B 
force and for Volkov No.2(m = 0.60) are very similar to each other in the region d > 2.5 fm. In the small-distance 
region d < 2 fm, the a-a potential in the case of the FSB force is shallower than that for Volkov No. 2 because of 
the three-body repulsive term. The a-a potential for the Gogny Dl force and the MVl case(l) force(m = 0.60) are 
similar to each other. Both forces indicate a smaller inter-cluster potential compared with the Volkov No. 2 (m = 0.60) 
force. This result is consistent with the a-a scattering behavior. 



B. I'^O-a system 



The total and potential energy of the ^^0-a system is shown in Fig. 3 and Fig. 4. The system is represented 

by the Brink-type wave function. Namely, and a described by the harmonic oscillator wave functions are located 
at {x,y,z) — (0,0, (i/5) and {x,y,z) — (0, 0,— 4rf/5), where d is the inter-cluster distance. The intrinsic energy of 
clusters is subtracted. The width parameter of and a is chosen to be the optimum width which gives the minimum 
energy of ^°Ne for each forces. As shown in Fig. 3, the potential energy is similar among three forces(Gogny, MVl 
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FIG. 2. Energy and potential between ''He-'*He. The potential energy(total energy) of the ''He-'*He system as functions of 
the relative distance (d) is shown in the lower(upper) panel. The effective forces are Gogny Dl, Volkov(m = 0.60), MVl cases 
l(m = 0.60) and FSB. 



TABLE II. The binding energy of the band head states(Oj and 1^) of the parity doublets in ^"Ne calculated with the 
AMD+GCM method. The adopted Hamiltonian is same as Eq.8. The FSB case(l) and the MVS case l(m = 0.60) are used as 
the central forces. The width parameter is chosen to be = 0.20 and v = 0.165, which minimise the binding energy of ^°Ne, 
for the FSB case(l) and MVS case l(m = 0.60) forces, respectively. Excitation energy of the Ij" states are also shown. 



(MeV) 

E.(ir) 



FSB 
157.2 
151.4 
5.8 



MVS 
15S.1 
146.5 
6.5 



exp. 
160.647 
154.859 

5.788 
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FIG. 3. Energy and potential between ^^O-'^He calculated with the Brink model wave functions. The potential energy(total 
energy) of the ^^O-^fie system as functions of the relative distance d is shown in the lower(upper) panel. The effective forces 
are Gogny Dl, MVl case l(m = 0.60) and FSB forces. Width parameters are chosen to be optimum widths which give the 
minimum energy for ^°Ne, i.e. v = 0.16,z/ = 0.165 and u = 0.20 for the Gogny Dl, MVl casel(m = 0.60) and F3B forces, 
respectively. 
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FIG. 4. The total energy of the ^®0-"'He system for the parity-projected states calculated with the Brink model wave 
functions, d is the relative distance between ^^O and *He. Width parameters are chosen to be optimum widths which give the 
minimum energy for ^°Ne, i.e. u = 0.165 and v = 0.20 for the MVl casel(m = 0.60) and FSB forces, respectively. 
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case l(m = 0.60) and F3B). In the total energy (Fig. 3), the minimum energy points exist at almost the same distance 
~ 4 fm for the MVl case l{m = 0.60) and FSB forces. In the small distance region, the FSB force has a more 
repulsive term of the total energy than the MVl case l(m = 0.60) force because the width parameter nu is larger 
for the FSB force than the MVl case l{m = 0.60) force, therefore, the kinetic energy increases more rapidly as the 
inter-cluster distance become small. 

The potential energy is closely related to the structure of parity doublet states in ^°Ne. The K = 0^ and 

the K = Oi bands in ^°Ne are known to be the parity doublet states, which are considered to originate from the 
cluster structure. In Fig. 4, the total energy of the parity-projected states for the Brink-type wave functions 
are plotted as a function of the inter-cluster distance d. Comparing the results with the FSB force and MVl force, the 
minimum energy is qualitatively similar. The minimum energy point for the positive-parity states shifts to the large 
distance region in case of the FSB force, which has an effect on the enhancement of the clustering in ^°Ne. 

In the above-mentioned discussion, the potential are studied within the framework, where the "'^^O and a are 

treated as inert clusters written by harmonic oscillator wave functions. However, in the real ^°Ne, the core excitation 
is more important than the case of a-a system. Moreover, the harmonic oscillator wave function is too simple to 
describe the real ^^O nucleus. In order to see agreements to the experimental data, we calculate the band-head 
statcs(0]'" and Ij^) of the parity doublets in ^"Ne with the AMD+GCM method, where the effect of core excitation 
is automatically included. The detailed description of the method is given later. As shown in Table II, the binding 
energy and the excitation energy of the Ij^ state are well reproduced by using the FSB force. 



C. Radii and Binding energy of Z = N(A < 40) nuclei and O isotopes 

In order to discuss the features of the saturation properties (radii and binding energy) in finite nuclei, we apply the 
FSB forces to microscopic calculations based on a AMD method for Z = N{A < 40) nuclei and O isotopes. z/ = 0.16, 
The formulation of AMD for a nuclear structure study of ground and excited states is explained in [1,2,17]. Here, we 
briefly explain the framework in the present work, namely, the generator coordinate method in the framework of AMD 
with constraint (AMD-I-GCM). The wave function of a system is written by superposition of AMD wave functions, 

$ = c^AMD + c'^'amd H • (5) 

An AMD wave function of a nucleus with mass number ^ is a Slater determinant of Gaussian wave packets: 

$AMn(Z) = ■^A{(pi,(p2, (fA}, (6) 
<i>x,(rj) oc exp[-iy{rj - ^)2], 



where the ith single-particle wave function (ipi) is a product of the spatial wave function (</>Xi), the intrinsic spin func- 
tion (xj J and the iso-spin function {xri)- The spatial part (0xj is represented by complex variational parameters(Xij, 
^2i, Xsj). Xii is the intrinsic spin function, defined by ^j, and Tj is an iso-spin function, which was fixed to be 
up(proton) or down(neutron) in the present calculations. Thus, an AMD wave function is expressed by a set of 
variational parameters, Z = {X„j, ^j} (n = 1, 2, S and i = 1, • • • , A), which stand for the centers of Gaussians of the 
spatial part and the intrinsic spin orientations of the single-particle wave functions. By using the frictional cooling 
method, we performed energy variation with respect to a parity-projected AMD wave function, while satisfying 
the condition that the principal oscillator quantum number of the system, {a'^a) = ' 
equals a given number. After variation with the constraint on the oscillator quanta, we superposed the spin-parity 
eigen states projected from the obtained AMD wave functions so as to diagonalize the Hamiltonian matrix and the 
norm matrix. The expectation values for such structure properties as the binding energy and the root-mean-square 
radius were evaluated by the superpositions. In the present calculations, we superposed several spin-parity projected 
AMD wave functions with respect to the generator coordinate, (a'^a). 
The used Hamiltonian in the AMD-I-GCM calculations was, 

^2 

H = ~ 2M 5Z ~ Vcoulomb + Vcentral + Vspin-orbit, (8) 

i 

where the first term stands for the kinetic energy operator, the second is for extraction of the center of mass energy, 
and Vcoulomb means the Coulomb energy term. The central force, Vcentral, is the FSB force. In addition to the central 
force, we employ the following two-body spin-orbit ioice{Vspin-orbit) '■ 
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TABLE III. Calculated binding energy and root-mean-square radii for a point-like proton distribution. Theoretical values 

for the FSB case (1) were obtained by AMD+GCM calculations. The values with the Volkov force No.2(V2), A modified 
Volkov force(MVl) case 1 and Tohsaki's Fl force are those calculated by the closed-shell HO wave functions ( [12,13]). The 
center-of-mass correction in the theoretical radii for V2, MVl, Tohsaki's Fl forces is subtracted as {r^) — 3/2aA {A is the mass 
number and a is the harmonic oscillator size parameter related to ?uj = h^a/M, where M is the nucleon mass.) The theoretical 
values of the Skyrme III force are three-dimensional HF-I-BCS calculations by N. Tajima et al [5]. The value in parenthesis was 
calculated with a single AMD wave function (the values for SIII are taken from [19]). The experimental r.m.s.r. were deduced 
from the charge radii [20]. 





Force 


V2(m = 0.65) 


MVl (easel) 


SIII 


Tohsaki's Fl 


F3B(this work) 


exp. 


*He 


E{MeY) 


28.0 


28.2 


26.5(25.6) 


27.5 


28.2(28.0) 


28.296 




r.m.s.r. (fm) 


1.46 


1.64 


1.88(1.66) 


1.50 


1.52(1.48) 


1.47 


160 


S(McV) 


105.9 


122.2 
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129.7 
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r.m.s.r. (fm) 


2.35 
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FIG. 5. Binding energy oi Z = N nuclei from 
with the experimental data. The theoretical values by HF-I-BCS calculations with SIII [5] are also displayed 



^spin— orbit — ^2<j \ (^) 

yiLS) ^ 1±P.1±P^,^. . . ^2^^ [_(-£.).] } , (10) 

(LS) (LS) 

where the same range of parameters as those of the G3RS force [18], a\ — 0.4472 fm and 03 = 0.5774 fm, are 

(LS) (LS) 

used. The strengths of the spin-orbit force were chosen to he v\ — —v^ = 1800 MeV in the present calculations. 

In Tabic III, the results of root-mcan-squarc-radii(r.m.s.r) for point-like proton distributions of double-closed nuclei, 
""■He, and "'*'Ca are listed. It is fomid that, overall, the FSB force well reproduces the radii and energy of these 
nuclei as well as Tohsaki's Fl force. It has been already known that it is impossible to systematically describe the 
nuclear saturation properties by Volkov forces. For instance, in the case of Volkov No.2(V2) with the Majorana 
exchange parameter m = 0.65, the radii of ^^O and ^"Ca are largely underestimated. On the other hand, the MVl 
force and the SIII force fail to qualitatively reproduce the compact size of a real ^He. A HF-I-BCS calculation with 
SIII gives an extremely large radius of an a particle because the center-of-mass motion is not exactly extracted in the 
model space. By an AMD calculation, where the center-of-mass motion was microscopically extracted, the calculated 
a radius with the SIII force was found to be 1.66 fm, which is as large as those with the MVl force and the SIII 
force. Therefore, we can conclude that these forces with zero-range density-dependent terms always overestimate the 
correct size of an a particle. 

In Figs. 5 and 6, the present results of the binding energy and the proton r.m.s.r. oi Z = N nuclei up to '*'^Ca with 
the FSB forces are compared with the experimental data. The present results well agree with the overall data in this 
nuclear region. 

Although the parameter sets, case (1) and case (2), of the FSB force give almost the same results iov Z = N 
nuclei, the results tov Z ^ N nuclei are different because the ratios of the inter-nucleon interaction for the isospin 
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FIG. 6. Root-mean-square radii of the proton density for Z = N nuclei from ''He up to ''"Ca calculated by AMD+GCM 
with the FSB case (1). The experimental data were deduced from the charge r.m.s.r. [20]. The theoretical values by HF-I-BCS 
with SIII [5] are also shown. 
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FIG. 7. Binding energy of neutron-rich O isotopes calculated by AMD-I-GCM with FSB case 1 and case 2. The HF-I-BCS 
calculations with SIII [5] and the experimental data are also displayed. 



T = 1 system to that for the T = system are differcmt between cases (1) and (2). The theoretical binding energy of 
neutron-rich O isotopes are compared with the experimental data in Fig. 7. The calculated results with both sets of 
interaction parameters agree reasonably with the experimental data. ^^O is known to be a particle-unstable nucleus. 
Concerning the systematics of the two-neutron separation energy, the neutron drip Une in O isotopes is reproduced 
by the case (2) force. Although the drip line is not described by the case (1) force, we can not conclude that the 
parameter case (1) is inappropriate at the present stage because the present model space seems to be insufficient to 
describe the detailed structure of ^'*0, as conjectured from a disagreement between the theoretical radius and the 
experimental one. 

The present results concerning the r.m.s.r of O isotopes quantitatively agree with the experimental data deduced 
from the interaction cross-section data [21], except for ^^O (Fig. 8). The reason for the failure in reproducing the 
enhancement of the radius near to the drip line is conjectured that the present model space is not sufficient to describe 
the detailed behaviour of weakly bound neutrons. 



D. nuclear matter property 

The matter properties were evaluated by plane waves in the Fermi sea. The equation of state(EOS) of the symmetric 
nuclear matter is illustrated in Fig. 9. The properties at the saturation point arc listed in Table IV. The SIII force, 
Gogny Dl force, Tohsaki's Fl force, FSB force give similar results for the saturation energy Eg/ A and the saturation 
density ps{ps = 2fcj/37r^, where fc/ is Felmi momentum at the saturation point). On the other hand, the MVl force 
can not satisfactorily explain the following empirical values: Eg/ A ^ —16 MeV and p., ^ 0.16 fm^'^(fc| ^ 1.33 fm^^). 

The small incompressibility of matter(iCoo)) namely, the soft-type equation of state(EOS) is given by such effective 
forces as the Gogny force and the SLy4 forces(a Skyrme parametrization) [11], while the FSB, SIII, Tohsaki's Fl, 
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FIG. 8. Root-mean-square(r.m.s.) matter radii calculated by AMD+GCM with FSB case 2. The experimental data are 
those deduced from the reaction cross sections taken from Ref. [21]. The theoretical values by HF+BCS with Sill [5] are also 
shown. 



MVl forces give large values of incompressibility, which correspond to the hard EOS (see table IV). Comparing EOS 
for the FSB force with that for the Gogny Dl force, it is fount that the EOS is similar to each other in the low-density 
region, while the matter energy for the FSB force rapidly increases in the high-density region. The incompressibility of 
nuclear matter has been discussed in relation to the energy of giant monopole resonances (GMR) [22-25]. Although a 
soft-type equation of state(EOS), like = 200 ~ 240, is preferred to quantitatively describe the GMR energy in HF 
analysis with Gogny and Skyrm-type effective forces [23,22], there still remains an inconsistency with the theoretical 
values suggested by a relativistic-mean- field analysis [25], where larger incompressibilities are predicted. On the other 
hand, the hard-type Tohsaki's Fl force reasonably reproduces the GMR energy of ^"Ca [IS]. 

It has been known that matter saturation can not be described by effective two-body central forces without any 
repulsive density-dependent or three-body terms, as can be seen in the EOS given by the Volkov force (Fig. 9). 
In other words, the density-dependent or three-body terms are essential to explain the nuclear matter saturation. 
In order to exhibit the mechanism of nuclear saturation, we decomposed the total potential energy into different 
contributions, E^^'^\ in each subspace (ST), defined by the projectors pf^"^) = {1 ± Pa)il ± Pr)/^ in the same way 
as is the analysis described in Ref. [9]. The decomposed potential energy of the FSB forces is shown in Fig. 10. The 
dominant parts of the potential energy at a saturation density of A:/ ~ l.S fm~^ are the singlet even(^£') and triplet 
even(^£) terms, while the singlet odd (^O) and triplet odd (^O) parts are minor. We notice that the part plays 
an important role in the saturation mechanism. This result is consistent with a tendency which was indicated in 
the nuclear matter calculations by Brueckner-Goldstone Theory using realistic nucleon-nucleon interactions [26]. It 
is very reasonable because the matter saturation is considered to come from a suppression of the tensor part of the 
bare nucleon-nucleon interaction in nuclear matter, because the Pauli principle acts to reduce the attractive tensor 
force due to a blocking of the intermediate states. Therefore, tensor suppression in nuclear matter is regarded to be 
phenomenologically simulated by the repulsive density-dependent terms in the effective central forces. In that sense, 
the repulsive density-dependent term should be dominated by the triplet even component, just as in the case of the 
tensor force. As a result, the triplet even component indicates the saturation behavior with an increase of density. 
Strictly speaking, also, the "'O part of the potential energy has an enhancement in the region kj > 1.5 fm"'^, because 
the three-body repulsive term is not pure ^E, but it also contains the ^O component because of the simple form of 
the three-body operator in FSB forces. Such an enhancement in a rather high-density region may have at least no 
significant effect on the structure of finite nuclei. 



IV. DISCUSSION 



A. Significance of finite-range three-body term 



The finite-range three-body terms in the phenomenological effective forces are important to describe the overall 
properties for nuclei for a broad range of masses. As mentioned in III, and also already argued in the pioneering work 



10 



TABLE IV. Saturation Properties for nuclear matter. The saturation energy, Felmi momentum, and incompressibility at 
the saturation point for various effective forces. 



Force 


MVl(casel m=0.60) 


SIII 


SLy4 


Gogny(Dl) 


Tohsaki's Fl 


F3B(this work) 


E/A{MeY) 
fe/(fm-i) 


-21.01 


-15.87 


-15.97 


-16.32 


-17.0 


-17.9 


1.48 


1.29 


1.33 


1.35 


1.27 


1.32 


ifoo (MeV) 


299 


356 


230 


228 


309 


390 
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FIG. 9. Equation of state for symmetric nuclear matter. The energy per nucloon is plotted as a function of the density. The 
effective forces are Gogny Dl, FSB, Volkov No.2(m = 0.65), and MVl case l(m = 0.60). 



by Tohsaki [13], the finite-range three-body terms are empirically essential to simultaneously reproduce the properties 
of Q-particlc(sizc and a-a phase shift) and nuclear matter saturation. Moreover, the finite-range feature is reasonable 
considering the original range of the tensor force, whose suppression effect in matter is simulated by the three-body 
terms. In order to explain the indispensability of the finite-range nature in three-body terms, we demonstrate the fatal 
failure of general effective forces with zero-range density-dependent terms in describing the correct size and energy of 
an a particle and the saturation density and energy of nuclear matter. We assume an effective central force, which 
consists of a two-body term with a two-range Gaussian form and a zero-range density-dependent term, as follows: 

2 2 

Vij = {I - m - mPaPr) i;iexp(-^) +W2exp(-^) + v^p" {nj) ■ 5{rij), (11) 
where cr = | or 1. (12) 

Coulomb force is omitted. By assuming the (Os)"* HO wave function for an a particle and plane waves in the Fermi 
sea for nuclear matter, we can easily obtain the optimum size(ra) and energy (£^0) of an a particle and the matter 
saturation density(ps) and energy(£s) as a function of the parameter set (t;i, ai, W2, 02, Vz,m) as rQ(vi, ai, z;2, 02, W3, m) 
, Ea{v\,ai,V2, a2,vs,m) , ps{vi,ai, V2, a2,vs,m) , and Es{vi,ai,V2,a2, V3, m)/A. We try to search for a parameter set 
which can reproduce the empirical data of these four values. Firstly, we put the following conditions: that the size and 
energy of an a particle arc equal to the experimental data, = 1.48 fm and = 29.253 McV. Here, the experimental 
input of the energy is a modified one with a Coulomb-force correction. Keeping these two constraints, we vary the 
parameters {vi,a\,V2,a2,V3,m) in the following ranges: 0.6 < oi < 2.5 fm, 0.5 < a2 < a\ fm, 100 < V3 < 4000 MeV 
fm^ (cr = 1/3) or 100 < 773 < 4000 McV fm^((7 = 1), 0.01 < m < 0.99, and investigate the behaviour of the matter 
saturation point {ps,Es/A), given by the various sets of paramters. We also put a physical restriction, vi < 0, that 
indicates that the long-range term must be attractive. The saturation points {ps,Es/A) are displayed in Fig. 11. The 
upper panel indicates the saturation points in the case a = 1/3, and the lower panel corresponds to the case cr = 1. 
Under the conditions of the a-particle properties, the saturation points can not reach the empirical saturation point 
{ps ^ 0.16 fm"*^, Es/A ~ — 16 McV) in both cases, because the allowed area of the saturation points does not include 
the empirical one. This means that there is no reasonable parameter set which is able to simultaneously describe the 
saturation properties of an a particle and nuclear matter. In other words, with "zero-range" density-dependent-type 
effective central forces, which explain the matter saturation point, it is difficult to describe the compact size of a 
real a particle. From these results, we can conclude that the finite-range terms are indispensable to simultaneously 
reproduce the properties of a-particle and matter saturation. 
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FIG. 10. Partial contributions E^^'^^ of the potential energy in each subspace of singlet even(^E), triplet 



'E), singlet 



odd(^O) and triplet odd(^O) defined by the projectors P^'^^^ = (1 ± Pa)0- ± •Pr)/4, where (ST)=(0,1), (1,0), (0,0) and 
(1,1), respectively. The calculated values of the partial potential energy per nucleon are plotted as functions of the Fermi 
momentum(fc/). 
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Saturation density Pg (fm-^ ) 

FIG. 11. Dots indicate the saturation energy and density in symmetric nuclear matter obtained by effective central forces 
with finite-range two-body terms and a zero-range density-dependent term, which satisfy the experimental data of the size 
and energy of an a particle. The square points correspond to the empirical saturation energy and density {ps = 0.16 fm~^, 
Es/A = —16 MeV). The details are explained in the text. 



B. Relation with s-wave nucleon-nucleon scattering 

It has already been known that bare nucleon-nucleon interactions have leading tensor terms and hard cores, both of 
which are not easily treated in the microscopic calculations of nuclear many-body systems. Instead, one usually adopts 
phenomenological effective nucleon-nucleon forces dominated by central forces with no hard cores. It is regarded that 
the tensor terms and hard cores are renormalized in the effective central forces. As already mentioned, the necessity 
of repulsive density-dependent terms of effective forces for describing the saturation properties originates from a 
suppression of the attractive tensor force in nuclear matter [26]. In other words, the density dependence of the tensor 
force suppression is phenomenologically simulated by density-dependent repulsive terms in the effective central forces. 
From this point of view, we demonstrate a link between the effective forces and nucleon-nucleon(N-N) scattering and 
discuss the role of the density-dependent repulsive terms. Firstly, we reproduce the s-wave N-N scattering phase shift 
by the central two-body forces in the FSB force. Second, we discuss the role of the residual part with the three-body 
term, which should be connected with the matter effect on the N-N interaction. 

We decompose the central part of the present effective forcc(F3B) into two terms. One is a two-body term, y^-body^ 
which can reproduce the s-wave nucleon-nucleon scattering behaviour, and the other is the residual term. Ay, which 
contains the three-body term, defined as follows: 

Kent™; = E.<, V;f + E.<,<. = y^-body + AV, (13) 

^'-''°''' = E.<,^.f +E.<,^.;-, (14) 

AF=-E.<,V'/,+E.<,<.^.S' (15) 

(16) 

where v' = 32.5 MeV. Other interaction parameters in V^^' and V^^j! are same as those of F3B force. The first term 
y2-body -^^ ^-j^^ two-body central force. The s-wave phase shift for the low-energy N-N scattering data is 

well described by this two-body term y2-body^ g^own in Fig. 12. Moreover, the calculated scattering lengths of 
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FIG. 12. ,s-wavo phase shift of nucleon-nucleon scattering calculated with the effective two-body central force ^""^^ defined 
in equations 13-16 of the text. 
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FIG. 13. Partial contributions iS*^^' of the nuclear matter potential energy for AV in each subspace defined by the projectors 
p(ST) ^ (i±p^)(i±p.^)/4. 



the singlet and triplet .s-wavc for .^^.^ _ 5.37(fm) and at = — 19.6(fm), which reasonably agree with the 

experimental values, = 5.39 and at — —23.7. Therefore, we consider this y^-body effective force between 

bare nucleons, where the tensor force in the realistic interaction is renormalized in the central term. Thus, the second 
term(Ay) in Eq. 13 can be understood in terms of the matter effect in finite nuclei and nuclear matter. Since the 
leading matter effect originates from a suppression of the tensor force, AV should be repulsive, and behaves as triplet 
even. In Fig. 13, we represent the partial contributions, E^^'^\ of the potential energy given by AV in each subspace 
(ST) of nuclear matter. At a saturation density of kf ~ 1.3 fm~^, the dominant component of AV is the triplet even 
part. This means that the density-dependent repulsive term(AF) actually takes the role of tensor force suppression 
in nuclear matter. 

As mentioned above, the central effective force(F3B) is decomposed into two parts, v^2-body ^^^^ principle, 
it is easy to understand the decomposition if y^-body jg g^me as the two-body term J2i<j of ¥2)Q force, and 
AV consists of only the three-body term of FSB. However, in the present analysis, y^-body jg ^.j^g modified one from 
the original V^^ , and AV^ involves a two-body term in addition to the dominant three-body term. We consider 

the reason as following way. The origin of the finite-range three-body term may not be the pure three-body effect 
but it also contains some component from the two-body central force. It is because of the restriction of the simple 
finite-range three-body operator, where we adopt a common range for the distances of all nucleon pairs in three-body 
term and omit spin-dependent terms for simplicity to apply to practical calculations with AMD method. In order to 
overcome the restriction, we consider that AV", which consists of the two-body and thee-body terms, indicates the 
matter effect and means the effective density-dependent repulsive force. 
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V. SUMMARY 



Wc proposed new effective inter-nucleon forces with a finite-range three-body operator for microscopic calculations 
of nuclear structure. It was found that the proposed forces are suitable to describe saturation properties over a 
wide mass number region from the a particle to nuclear matter. Compared with similar-type forces proposed by 
Tohsaki, an advantage of the new forces is that they are more practical because of the simple three-body operator, 
and therefore are applicable to microscopic calculations of general nuclei. The proposed effective forces were applied 
to AMD-I-GCM calculations of Z = N {A < 40) nuclei and O isotopes. In addition to the success in systematically 
reproducing the radii and energy of these nuclei, we demonstrated the characteristics of the new forces concerning 
the a-a scattering behaviour and matter EOS. Especially, the new forces successfully describe the compact size of a 
real a particle and the a-a scattering phase shift, which typical effective forces with zero-range density-dependent 
terms can not reproduce. Concerning the matter properties, we found that the saturation mechanism in the nuclear 
matter originates from the dominant triplet even component in the repulsive three-body term. The new forces were 
also discussed in relation with the s-wavc N-N scattering behavior. 

The density-dependent repulsive terms are essential to describe nuclear saturation. The origin is a suppression of 
the tensor force of bare nucleon-nucleon interactions in nuclear matter. From this point of view, we discussed the role 
of the finite-range three-body term in the new forces. The importance of the finite-range three-body term was also 
described based on empirical concepts. We conclude that the finite-range three-body operator plays an essential role 
in describing the overall properties over a wide mass number region as well as the matter saturation properties. 
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